X 



Band structure and Bloch states in birefringent ID 
magnetophotonic crystals: An analytical approach 



Miguel Levy, Amir A. Jalali 

I Michigan Technological University, Department of Physics, 

^ ; 1400 Townsend Dr, Houghton, MI 49931, USA 

■ An analytical formulation for the band structure and Bloch modes in elliptically 

\^ , birefringent magnetophotonic crystals is presented. The model incorporates both 

(N 

the effects of gyrotropy and linear birefringence generally present in magneto-optic 

Q ' thin film devices. Full analytical expressions are obtained for the dispersion relation 

•i-H ' 

-4— > < 

pi |. and Bloch modes in a layered stack photonic crystal and their properties are 

o , 

I analyzed. It is shown that other models recently discussed in the literature are 

O 

* ^ ' contained as special limiting cases of the formulation presented herein. 

^ I. INTRODUCTION 

o 

, Photonic crystals in magnetic systems have been t 



ne suD 



ect o: 



a number of theoreti- 
sys- 



O ' cal and experimental studies in recent years. Q, I2I, 1^, ^, 0, 6, Q, [s] Interest in these 
O . tems stems from the scientific and technological possibilities resulting from the combina- 
O ! tion of magnetic tunability, non-reciprocal effects such as Faraday rotation, and band gap 



7|, 81 Examples of interesting phenomena that have been investigated 



engmeermg. 

p |. include electromagnetic unidirectionality, frozen light generation, the magnetic tuning of 



photonic bandgaps, and Faraday rotation 



enhancement 5, 0, fl H, O, Q, 11] 



Work reported in the literature on such photonic crystals has dealt with systems in which 



at least one of the constitutive component materials possesses circular birefringence. [12 1 
These are called gyrotropic or bigyrotropic materials. However, technological imperatives 
such as miniaturization, the development of integrated devices and optical waveguide fabri- 
cation demand that polarization effects resulting from stresses and optical confinement be 
taken into account. Magneto-optic films grown by various deposition techniques and sub- 
ject to lattice mismatch and differential thermal expansion with the substrate will acquire 



linear birefringence. [3, Q, [l^ In addition, in optical channels formed on the film, transverse 
electric (TE) and transverse magnetic (TM) modes will generally possess different effective 
or modal refractive indices, hence linear birefringence. 15(] The combination of circular and 
linear birefringence results in elliptically polarized modes, as pointed out, for example, in a 



recent publication by one of us. 16| Hence we prefer to use the term elliptical birefringence 
in this case. 

These phenomena have an impact on the polarization state of light propagating in a 
magneto-optic photonic crystal and on its band structure. Although some work has been 
done on the transmittance and polarization response in magnetic photonic crystal waveg- 
uides, the effect of elliptical birefringence on the band structure and Bloch modes in such 



systems has not been directly studied. [8|, llTI] Recently Khanikaev and co-workers presented 
an analytic solution to the band structure of a photonic crystal model consisting of ultra- 
thin magneto-optic layers in a stack Circular birefringence was included in the model 
but not the elliptical birefringence that results from combined optical gyrotropy and stresses 
or optical confinement. In the present article we present an analytic solution for the case 
of a gyrotropic photonic crystal stack with elliptical birefringence. The model does not as- 
sume ultra-thin magneto-optic layers and contains the solution obtained by Khanikaev and 
co-workers as a special case in the limit of zero linear birefringence and ultra-thin magnetic 
layers. In addition, the model does not confine itself to quarter-wave plates but allows the 
layers to have arbitrary thickness limited only by the period of the crystal. The model 
supports the combination of magnetic and non-magnetic layers as well. 

The value of an analytical solution to the case of birefringent magnetophotonic crystals 
is that it gives explicit form to the dependence of the band structure and the Bloch states 
on the material and geometrical parameters of the crystal. It thus represents a useful tool 
in the analysis and design of photonic band gap systems. The formulation presented here 
shows that the wavevector dependence of the optical frequencies is sensitive to the relative 
birefringence between adjacent layers in the photonic crystal stack. And it gives explicit 
functional form to this dependence for Bloch states of different helicities. Layer thickness 
and thickness ratios between different layers also strongly impact the band structure, as 
shown in the discussion on Bloch modes and the dispersion relation. 
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II. WAVE PROPAGATION IN BIREFRINGENT MAGNETOPHOTONIC 

MEDIA 



The wave equation in a birefringent magneto-optic material is given by 181] 

(fc^e- fc^I + kk) ■ Eo = 0. 



(1) 



where fco = cu/c, c is the speed of hght in vacuum, Eq is the plane wave amplitude, kk is 
a dyadic product of the wave vector, and I is the identity matrix. At optical wavelengths, 
the relative permeability fi differs only slightly from unity whereas the relative permittivity 
tensor e has the form 
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(2) 



Here the saturation magnetization of the magnetic medium is directed along the z-axis and 
rotational symmetry about this axis is assumed. Linear birefringence is expressed through 
the difference in diagonal components e^x and e^^, and gyrotropy through the off-diagonal 
elements iexy and —iexy Our treatment assumes a monochromatic time dependence exp(ia;t) 
and no optical absorption. By solving the wave equation in this birefringent magneto-optic 
medium, one obtains the eigenvectors 

^ cos a ± sin a ^ 

1 
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cos a ± sin a 
±i cos a — i sin a 




(3) 



with refractive indices n±, and = e ± + e^y, respectively. Here e = {eyy + exx)/'2, 
A = {eyy — e^,^.)/2, and a, referred hereafter as elliptical birefringence parameter, defined as 



a 



7/2, with 7 given by 
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In deriving Eq. ([3]), wave propagation parallel to the 2;-axis is assumed (Fig. [T]). Expres- 
sions ([3]) for the eigenvectors are identical to the elliptical eigenmodes introduced in Ref. [l^ 
for elliptically birefringent magneto-optic media. These eigenvectors can be simply related 
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to circularly polarized modes through a rotation transformation on c_|_ and c 

(5) 




cos t: sin 77 



sin ^ cos 2 




where 



III. ONE-DIMENSIONAL BIREFRINGENT MAGNETOPHOTONIC CRYSTAL 

The basic geometry of our model is depicted in Fig. ([1]). Plane waves are normally incident 
on a periodic stack structure consisting of alternating elliptically birefringent magneto-optic 
layers. The layers have different average dielectric constants e, and are not assumed to 
have the same linear birefringence terms A, or gyrotropic components e^y In addition, the 
model does not assume quarter-wave plates but allows the layers to have arbitrary thickness. 
Thus the model is quite general and does not impose any constraints on the relative linear 
birefringence, gyrotropy, or thickness of the layers. The electric field in the n-th layer of the 
periodic stack is a linear combination of the elliptical eigenvectors obtained in the previous 
section. 



Eoi exp (^-i'^n+{z - Zn)j + ^02 exp (^i'^n+{z - z„) j e+ (7) 
-EosGxp (^—i—n_{z — Zn)^ + £'04 exp (^i—n_{z — z^ 



c / \ c 

From Eq. we notice that it is essential to perform the transformation introduced in 
Eq. (I5l) to circularly polarized representation, since the eigenvectors e± change across the 
boundary. This concerns the continuity requirement of the tangential components of electric 
and magnetic fields at each interface. Following Yeh's matrix formulation [l^, the transfer 
matrix which relates the four eigenmode amplitudes in the (n — l)-th layer to the amplitudes 
in the n-th layer is obtained as (Fig. [T]) 
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where the dynamical and propagation matrices D and P are defined as 
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sma 



,(n) 


cos a 


— sin a'-"^ 


— sin a*^"-* 




— cosa^ -^n^. 


(n) (") 

— sma^ m_ 


(n) (") 

sm a*- -'n_ 




sin a'-"-' 


cos a'-"'' 


cosa*^"^ 




(n) {") 

— sm 'n . 


(n) (") 

COS >n_ 


(n) (") 

— cosa'^ 



(9) 



p(n) 








(n) 



with 





e" 

\ e" ^ 

where d^"-* is the thickness of the n-th layer. 



\ 






(10) 



From Eqs. (!8l)- (fT0|) we note that for zero birefringence the transfer matrix becomes block 
diagonal and reproduces the results obtained in Ref. 201] for circularly polarized eigenmodes. 



IV. BLOCK MODES AND DISPERSION RELATION 

Consider a unit cell of the birefringent magnetophotonic crystal. This consists of two 
adjacent layers as shown in Fig. [H According to the Floquet-Bloch theorem, the optical 
electric field obeys the relation 



where Ek{z) is a periodic function of coordinate z with a period A. 
From Eq. ([8]) the transfer matrix for a unit cell of the crystal obeys 

^ V ' 

S(") 



The Floquet-Bloch theorem thus imposes the following condition: 

rp(n-l,n+l)ji ^ ^j.^ 



2l| 



(11) 



(12) 



(13) 



To solve this eigenvalue problem we note that the term S*^"-* in Eq. ( |T2|) can be transformed 
into block diagonal form through a similarity matrix transformation as follows 



g(n) _ (XjW) ^S(")U("). 



(14) 
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where U^"-* and Si"^ are given by 



Ijin) 
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respectively. Upon substitution of Eq. ffT^ into Eq. f|T2l) . we obtain 
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rj-i(n- 1,71+1) _ ^^(n,n+l)^ g(n)^(n,n+l)p("+l) 
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The matrix is given by 
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(15) 



(16) 



(17) 



(18) 



, (19) 



where x^"'""*"^^ = — a*-"'^^^ We notice that only the relative elliptical birefringence 
parameter x*-"'"'^^-' appears in the unit cell transformation matrix. 

The advantage of reformulating the problem in this way is that it allows us to express the 
Floquet-Bloch condition in terms of block diagonal matrices, making it possible to obtain 
fully analytic solutions. Upon substitution of Eq. f|T71) into Eq. f|T3l) we obtain a generalized 
eigenvalue equation of the form 



(20) 
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where 8 = $("-'"-+i)p"+^E. The eigenvalues of the transformation matrix T*^" are 
given by the solutions to 



ISl") - ABl 



0. 



The elements of the matrix B are given by 
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(21) 

(22) 
(23) 

(24) 

(25) 

(26) 
(27) 
(28) 

(29) 
(30) 



Notice that all off-block diagonal components contain a difference term between trigono- 
metric functions of (3+ and /?_. For all available magnetic garnet materials in the infrared 
wavelength range, and most other Faraday rotation materials of interest such as Fe:InGaAsP 



and CdMnHgTe 



22 



23| . this difference is very small. For exarnple, a typical value for the 



dielectric tensor of bismuth-substituted lutetium iron garnet is IGj 



v 



5.369 i0.00274 
-i0.00274 5.373 

e,. 



(31) 



yielding, for one-quarter and one-third wave plates, the following B matrices 

/ -0.000681 -i0.431 -0.000196 i0.0000539 \ 

-i2.32 -0.000681 -i0.00029 -0.000196 

-0.000196 i0.0000539 0.000681 -iO.432 

\-i0.00029 -0.000196 -i2.32 0.000681 / 



(32) 
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(33) 



and 

/ -0.501 -zO.373 -0.000227 z0.000103 \ 
-z2.01 -0.501 i0.0000525 -0.000227 
-0.000227 i0.000103 -0.499 -iO.374 
\i0.0000525 -0.000227 -i2.01 -0.499 / 

respectively. In both cases we can see the off-block diagonal elements to be much smaller 
than the block diagonal elements. This is also true for technologically important non- 
garnet materials such as Fe:InGaAsP with Verdet coefficient of 12.5° / mm T and e = 9.61 
at a wavelength of 1.55 fim, where similar order magnitude values for the B matrix are 
obtained. 22| In Fig. [2] we show the variation of four elements of the B matrix, two elements 
from the block diagonal and two off-block diagonal elements, versus layer thickness d^^'' 
assuming the dielectric tensor Eq. fl3ip . For any value of layer thickness {S^^ < wavelength) 
the off-block diagonal terms are small. 

Moreover, all off-block diagonal elements, themselves already small, will be raised to the 
second or fourth powers in the determinant of the 4x4 matrix making their contributions 
even smaller. Therefore, we can safely ignore the effect of off-block diagonal terms on the 
determinant and set them equal to zero. The problem now reduces to an eigenvalue problem 
for two 2x2 matrices. The solution for the A's yields the Bloch phase factor exp{iKA) and 
the dispersion relation, 

cos KA = B± cos - ^C± sin /^i"^ , (34) 
with B^{BJ) = i?i 1(533). The factors C± in Eq. fl34|) are defined as 

C± = zni"^5,,, + 4^5,,., (35) 

with the subscripts {i,j} take {1,2}({3,4}) for {+}({ — }) sign. Equation flMl) determines 
the dispersion relation for the Bloch wave vector K with frequency. 

If the linear birefringence term A in this expression is set equal to zero, and the unit cell 
is composed of a thin magneto-optic layer and an isotropic layer, Eq. flMl) reduces to the 
ultra-thin approximation of Ref. ^ upon letting the thickness of the magneto-optic layer 
tend to zero. In this limiting case Eq. flMl) becomes: 



cos KA ^ cos koA — - ± 2 ^'^ k^A, (36) 



in agreement with Ref. j^. Here A = (ra^ — in?i) d and m = n^d are parameters used in the 
above referenced work. The present formulation also agrees with the treatment in Ref. [l^ 
that discusses the limit of uniform linear birefringence across the photonic crystal and where 
the difference between refractive indices of the two layers in the unit cell is small. In that 
case a'-"^ = a'^""'"^) and dispersion relation Eq. ( IMl) reduces to 



cos KK = cos (z?!"^ + A 



(n+l) 



(37) 
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with identical elliptical eigenmodes and average indices n± to those in Ref. IGj. 
The eigenvectors in Eq. (I13p are 
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V 1 / 

This model yields elliptically polarized solutions with different elipticity across layer bound- 
aries, according to Eq. ([7]). From the eigenmodes Eq. fl39l) we notice that the elliptically 
polarized electric field of a layer in a unit cell not only depends on the elliptical birefrin- 
gence parameter a*^"-* of that layer but also on the relative elliptical birefringence parameter 

^(n,?i+l) 

Band structures based on dispersion relation Eq. fl34l) for two crystals having different unit 
cells are plotted in Figs. ([3]) and (jl]). The unit cell consists of one birefringent [Bi,Lu]3Fe50i2 
layer adjacent to a Lanthanum garnet layer (La3Ga50i2) with the dielectric constant given 
24( 1 . Figure ([3]) corresponds to the case d*^") = d'^"''^^^ = 0.5A, and Fig. (jl]) to 



in Ref- 

= 0.3A. 

Although the gyrotropic components and linear birefringence term A are naturally 
small we are able to tell apart the two branches of the electromagnetic modes in the unit 
cell especially in the higher Bragg orders. Of particular note is the fact that the relative 
thickness of layers d^"''' and d^''^^^^ strongly affects the band structure, as can be seen by 
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comparing the third and fourth-order Bragg mode branches in Fig. [3] and Fig. HI Although 
the material compositions are the same in both figures, a departure from the quarter-wave 
plate configuration in Fig. [3] shifts the band gaps enough to induce a cross-over between 
different helicity modes in the high-order branches. We have also calculated numerically the 
band structures using the full form of the B matrix in eigenvalue Eq. (l20l) and compared with 
the analytic expression Eq. The maximum difference between the two formulations is 

at most 0.3%. 



V. CONCLUSIONS 



We have presented a fully analytical treatment of the band structure and Bloch states 
for one-dimensional elliptically birefringent magnetophotonic crystals. Expressions for the 
dispersion relation and Bloch modes are obtained. This formulation is applicable to most 
technologically important magneto-optic systems such as iron garnets in the near infrared 
and diluted magnetic semiconductors in high and low magnetic fields. The model treats the 
case of photonic crystals in magnetic systems in the presence of non-circular birefringence, 
such as result from strains due to lattice mismatch and differential thermal expansion with 
the growth substrate. In other words, the formulation presented here accounts for gyrotropy 
and linear birefringence simultaneously. This is particularly noteworthy because of the 
large class of photonic bandgap structures on magneto-optic films encompassed under this 
category and that naturally arise from standard fabrication techniques. 

A key feature of the treatment is the reduction of the Floquet-Bloch eigenvalue equation 
to block diagonal form, and hence its analytic solution. The model reduces to the exact 
solution for one- dimensional magnetophotonic crystals with ultrathin magnetic layers in the 
absence of linear birefringence. It also contains the solution for magnetophotonic crystals 
with uniform birefringence as a limiting case. Both of these models have been recently 
discussed in the literature. [3], IigI The power of the formulation presented here is that it 
applies to layers of arbitrary thickness and elliptical birefringence level and that it allows 
for alternating birefringence from one layer to the next. 

Band structures and the form of the Bloch wavefunctions are found to depend explicitly on 
the relative birefringence parameter . Different branches arise in the band structure 

separating Bloch modes with positive and negative helicity. It is seen that the frequency 
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shift between these branches is strongly dependent on the relative thickness of adjacent 
layers in the stack. The band gaps grow with Bragg mode number generating cross-over 
points between different helicity states in high order for asymmetric structures. 
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FIG. 1: Schematic diagram of a one-dimensional birefringent magnetophotonic crystal with period 
of A. The magnetophotonic crystal extend indefinitely in the x and y directions. A plane wave is 
incident from the left onto the structure. 



15 




FIG. 2: Variation of elements of the B matrix versus layer thickness d*^") in the unit cell for a 
typical [Bi,Lu]3Fe50i2 layer. Solid lines with circles and dots are curves for Bu and B21 from 
block-diagonal elements, respectively. Dashed and dotted lines are curves for B23 and B13 from 
off-block diagonal elements, respectively. 
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FIG. 3: Band structure for a periodic stack of Bi:LuIG and LaGG with d^") = 0.3A and = 
0.7A, respectively. The dashed Une shows the wave branch and the black line shows the K_ 
wave branch. 
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FIG. 4: Band structure for a periodic stack of Bi:LuIG and LaGG with = = 0.5A. The 

dashed Une shows the wave branch and the black hne shows the K_ wave branch. 



